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ABSTRACT 

Context. The sample of Hot Jupiters that have been studied in great detail is ever-growing. In particular, when the planet transits 
its host star, it is possible to measure the planetary radius and (with radial velocity data) the planet mass. For the study of planetary 
atmospheres, it is essential to obtain transit and occultation measurements at multiple wavelengths. 

Aims. We aim to characterize the transiting hot Jupiter WASP- 19b by deriving accurate and precise planetary parameters from a 
dedicated observing campaign of transits and occultations. 

Methods. We have obtained a total of 14 transit lightcurves in the r'-Gunn, IC, z'-Gunn and I+z' filters and 10 occultation lightcurves 
in z'-Gunn using EulerCam on the Euler-Swiss telescope and TRAPPIST. We have also obtained one lightcurve through the narrow- 
band NB 1 1 90 filter of HAWK-I on the VLT measuring an occultation at 1.19 \im. We have performed a global MCMC analysis of all 
new data together with some archive data in order to refine the planetary parameters and measure the occultation depths in z'- band 
and at 1.19 [xm. 

Results. We measure a planetary radius of Rp = 1.376 ± 0.046 Rj, a planetary mass of Mp = 1.165 ± 0.068 Mj, and find a very low 
eccentricity of e = 0.0077^q55q32, compatible with a circular orbit. We have detected the z'-band occultation at 3(T significance and 
measure it to be dFocc.-J = 352 ±116 ppm, more than a factor of 2 smaller than previously published. The occultation at 1.19 |,im is 
only marginally constrained at SF„cc.NBim = 1^1 PP™- 

Conclusions. We have shown that the detection of occultations in the visible is within reach even for Im class telescopes if a con- 
siderable number of individual events are observed. Our results suggest an oxygen-dominated atmosphere of WASP-19b, making the 
planet an interesting test case for oxygen-rich planets without temperature inversion. 
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1. Introduction 

At the time of wiiting, about 290 planets have been confirmed 
to be transiting in front of their parent star^ Via the precise 
measurement of transit lightcurves, we are able to constrain the 
planetary radius, orbital inclination, mass (usually with the help 
of radial velocity measurements) and hence the planetary den- 
sity. 

Transiting planets open up a window into the study of 
planetary atmospheres, their structure and composition. High- 
precision spectroscopic or spectro-photometric observations of 
planetary transits allow to search for wavelength dependencies 
in the effective planetary radius, and from there conclude on 
the molecular species present in the planetary atmosphere. Also, 
from the transit lightcurve an independent measurement of the 
stellar mean density can be obtained ( Seager & Mallen-Omelas 
2003| l, which is particularly useful since it can be used to refine 



* Based on photometric observations made with HAWK-I on the 
ESO VLTAJT4 (Prog. ID 084.C-0532), EulerCam on the Euler-Swiss 
telescope and the Belgian TRAPPIST telescope, as well as archive data 
from the Faulkes South Telescope, CORALIE on the Euler-Swiss tele- 
scope, HARPS on the ESO 3.6 m telescope (Prog. ID 084-C-0185) and 
HAWK-I (Prog. ID 083.C-0377(A)). 

** The photometric time series data in this work are only available 
in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) or via|http://cdsweb.u-strasbg.fr/cgi-bin/q cat?J/A+A/| 
' Based on www.exoplanet.eu |Schneider et al.[2oTT| 



the host stars parameters, namely its radius and mass. Having 
more accurate knowledge of the stellar parameters directly trans- 
lates into more accurate physical values for the planetary mass 
and radius. This has lead to several campaigns collecting transit 
lightcurves of pubUshed planets, as done by e.g. Holman et al. 
(2006) and 'South worth et al. ( 2009| . As the parameter measured 
from transit lightcurves is not the planetary radius itself but in 
fact the dimming of the star by the planetary disc, these mea- 
surements are affected by the brightness distribution along the 
stellar disc, i.e. stellar limb darkening as well as occulted and 
non-occulted spots. Occulted dark spots or bright faculae lead to 
short-term flux variations during the transit as the planet passes 
areas of the star having a different temperature and thus a dif- 
ferent brightness. Non-occulted spots alter the stellar brightness 
outside of the planet's path, leading to a slight change in the ob- 
served transit depth. Depending on the spot distribution on the 
stellar surface, spots cause a rotational modulation of the stellar 
flux, with typical amplitudes of a few percent in the optical and 
timescales of several days. While the effect on the transit depth 
is small (100 ppm for a typical brightness variation and transit 
depth of both 1 %), it is in the order of precision needed for the 
detection of elements through transmission spectroscopy. Next 
to these physical effects, ground-based photometric lightcurves 
are known to suffer from correlated noise due to airmass, seeing, 
or other external variations (Pont et al. 2006). These effects can 
be mitigated by choosing optimal observation strategies (such as 
staying on the same pixels during the whole observation and de- 
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focusing in order to improve the sampHng of the PSF) but can 
rarely be completely prevented. In this work, we collect a large 
number of transit lightcurves from Im class telescopes and com- 
bine them in order to find not only a very precise but also accu- 
rate measurement of the overall transit shape. 

Observing the occultation of a transiting planet 
( Charbonneau et al. 2005 Deming et al. 2005 1 allows us 



In this paper we present recent observations of WASP-19 ob- 
tained both in optical and IR light. Our dataset includes 14 tran- 
sits and 10 occupations observed at optical wavelengths using 
EulerCam (Lendl et al. 2012 ) and TRAPPIST ( |GiUonetal.|20lT| 



Jehin et al.,2011|) and one oc cultation obtained at 1.19 \im with 



to measure the brightness ratio of planet and star and hence 
measure the flux emitted or, at shorter wavelengths, reflected, by 
the planet. At optical wavelengths, occupations have been mea- 
sured almost exclusively from space ( [Alonso et al. 2009 Snellen 



HAWK-I ( IPirard et al.|2004l|Casali et ain2006| ). We describe all 
observations and their reduction in Section|2|and give details on 



the modeling in Section [3] In Sections |4] and |5 
discuss the results before concluding in Section 



we present and 



^ral.||2009| IBorucki et'al.||2009) using the CoRoT and Kepler 2. Observations and Data Reduction 



satellites. Ground based observations have been few (jSing & 



L6pez-Morales||2009l |L6pez-Morales et al^|2010| [Smith et al 
201 1) and so far none of the detections has been independently 
confirmed. Observations of occupations in the infra-red have 



been plentiful, both from space (starting with Charbonneau et al. 
2005 1 |Deming et aL 20051 and from the gr ound, e.g. |de Mooij 



& Snellen] ( |2009D , |Gillon et~aL] ( |2009l l and |CroU et al.| ( |2010| l. 



These observations provide information on the composition 
and the temperature profile of the planetary atmospheres. Some 
planets, such as HD209458 show a temperature inversion at 
high altitudes ( Knutson et al.|2008]), that is usuall y attributed to 
high abunda nces of TiO and VO ( [Hubeny et al.||2003| |Fortney| 
|et al.|[2008l l. These molecules are efficient absorbers of stellar 
radiation heating up the high altitude atmosphere. However 
it is not yet clear why some planets show inversions while 
others do not. As the number of planets with characterized 
atmospheres increases, the presence of inversions is turning 
out not to be dependent only on either the incident stellar flux 



( Fortney et al. 2008,) or the host star activity level (Knutson 
et al.|2010| l. Recently, Madhusudhan et al. ( 2011| l suggest an 



additional connection between the C/O ratio and the presence of 
an inversion, as in atmospheres dominated by carbon, the main 
absorbers TiO and VO are not abundant enough to cause an 
inversion. It is essential to increase the sample of well-studied 
transiting hot Jupiters, and to provide accurate values for the 
measured occultation depths. 

WASP- 19b has been identified as a hot Jupiter by 

(|2010| l, based on data taken by the WASP survey (jPa 

|2006| l' The Mp = 1.1 7My and slighdy bloated (p = 0.44 pj) 
planet is orbiting a magv - 12.3 G8V dwarf with a period of 
0.79 days. At this close separation, the planet is assumed to 
have been undergoing orbital decay moving it to its current or- 
bital position at 1.21 times the Roche limit ( Hebb et al.||2010 



Hebb et al. 



lacco et al. 



Hellier et al.||201 1) . The star is known to be active showing a 



rotational modulation with a period of 10.5 days in the discov 
ery lightcurves ( Hebb et al.||201 0). Also, anomalies in transit 
lightcurves attributed to spot crossings have been reported by 
[Tregloan -Reed et al.| ( |2012| l. The projected stellar rotation axis of 
WASP-19 is ahgned with the planet's orbP pellier et al.|2011| l. 

Occupations of WASP- 19b have been measured in the past 
by |Gibson et al7| ( j2010| and |Anderson et al.| ( [20T0l ) usin g HAWK- 



I in the K and H band, respectively, as well as by Anderson 
|et al.| ( [20TT] ) using the Spitzer Space Tele scope at 3.6, 4.5, 5.8 



and 8.0 (xm. Recently, IBurton et al. ( 2012| ) have published a z'- 
band lightcurve obtained with ULTRACAM during one occul- 
tation of WASP-19b, claiming its de tection at 880 + 119 ppm. 
From the ensemble of measurements, [Anderson et al.| ( |201 1 j ) and 
Madhusudhan ( 2012[ l determine that WASP- 19b does not posses 
a temperature inversion. Using the z'-band value of [Burton et al ' 



( 2012| l, models favor a C-rich atmosphere. For the eccentricity 
of WASP- 19b, [Anderson et al.| ( [20n i derive 3-cr upper limit of 
e < 0.027. 



Between May 2010 and April 2012, we obtained a total of 25 
lightcurves of WASP-19. 14 of these observations were timed 
to observe the primary transit, while 11 were performed dur- 
ing the occultation of the planet. We made use of three in- 
struments: EulerCam at the 1.2 m Euler-Swiss Telescope and 
the automated 0.6 m TRAPPIST Telescope at ESO La SiUa 
Observatory (Chile), as well as HAWK-I at the VLT/UT4 at ESO 
Paranal Observatory (Chile). We include in our analysis also the 
Faulkes South Telescope (FTS) lightcurve published by [Hebb 



let al.l (|2 010| l, the HAWK-I H-band observation by [Anderson 
et al. ( 20l0|l and the radial velocity measurements presented in 



[Hebb et al.| ( [2010| l and |Hellier et al.| ( [201 l| l. All new observations 
are summarized in Table [IJ 



2.1. EulerCam 

Five transit and six occultation lightcurves have been obtained 
with EulerCam, the imager of the Euler-Swiss telescope at La 
Silla. The instrument and the reduction of EulerCam data are 



described in detail by [Lendl et al. (2012j). The observations 
were done either with a focused telescope or applying a small 
< 0. 1 mm defocus yielding stellar PSFs with a typical full width 
at half-maximum (FWHM) between 1.1 and 2.5 arcsec, while 
the exposure times were between 60s and 120s, depending on 
filter and conditions. On 2011-03-12, the CCD temperature dur- 
ing the observations was slightly elevated, -100° C instead of 
the nominal -115° C. The lightcurves were extracted using rel- 
ative aperture photometry, with the reference stars and apertures 
selected independently for each observation. 

2.2. TRAPPIST 

Nine transits and four occupations of WASP- 19b were observed 
with flie robotic 60 cm TRAPPIST (Gillon et al. 2011, Jehm| 
et al. 201 1) that is also located at the La Silla site. We defocused 
the telescope slightly in order to spread the light over more pixels 
yielding typical FWHM values between 3.6 and 6.5 arcsec on the 
images and used exposure times between 15s and 40s. Again, 
the lightcurves were obtained with relative aperture photometry, 
where IRAFj^is used in the reduction process. 

2.3. HAWK-I 

We obtained photometry of WASP- 1 9 using the HAWK-I instru- 
ment dPirard et al.|[2004| [Casali et al.|[200ll l on VLT UT4 dur- 
ing two occupations of WASP- 19b using the narrow band filters 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1: All transits observed together ■with their models and residuals. Instrument and filter used are indicated in each frame. All data 
have been binned in t'wo minute intervals. 
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Table 1: A summary of newly obtained photometry. Date, instrument, filter and the nature of eclipse are given for each observation 
together with the photometric model function, red noise amplitude Pred (as defined in |Winn et al.| ( |20d8l )) and the RMS of the binned 
(5 minutes) residuals. The notation pij) refers to a polynomial of degree i of parameter j, e.g. p(f^) denotes a polynomial of second 
degree with respect to time. 



NBl 190 and NB2090. Unfortunately, during the NB2090 obser- 
vations, the target exceeded the linearity range of the detector, 
as a consequence the data do not have the necessary precision to 
detect the occultation. However the NBl 190 data are good, and 
thus we restrict our analysis to them. 

The NBl 190 observations of WASP-19 took place on 20 
January 2010 from 02:55 to 06:55 UT covering the predicted 
occultation time together with 145 min of observations outside 
of eclipse. The detector integration time (DIT) was kept short 
(3s) so the counts of the target and a bright reference star did not 
exceed the linear range of the detector. The data were obtained 
alternating between two Jitter positions, in order to be able to 
coiTect for backgi-ound variations if necessary. 

The data were corrected for dark and flat field effects us- 
ing standard procedures. Then, we identified bad pixels on the 
images and substituted their values by the mean of the neighbor- 
ing pixels. Here, we experimented with different cutoffs for the 
identification of bad pixels and obtained the best results discard- 
ing pixels deviating by 4 cr for background values and 40 cr for 
stars. The target flux was extracted from the corrected images 
using aperture photometry. We tested a set of constant apertures 
as well as apertures that varied from image to image as a func- 
tion of the stellar FWHM. The sky annulus was kept constant for 
all images. The best result was obtained using a variable aper- 
ture of two times the FWHM. We tested all bright stars of the 
four HAWK-I chips and found the best photometry using only 
the single bright star located on the same chip as WASP-19. The 
bright stars on the other detectors showed significantly different 
variation and were thus not used. The lightcurves are shown in 
Figure |2] 



3. Modeling 

We performed a combined analysis of all photometric (transit 
and occultation) data together with published radial velocities. 
The data were modeled using the Markov Chain Monte Carlo 
(MCMC) method in order to derive the posterior probability 
distributions of the parameters of interest (see 3.1 for details). 
Incorporated in our analysis are models for photometric correc- 
tion functions which account for photometric variations not re- 
lated to the eclipse, i.e. airmass, weather or instrumental effects. 
We also rescale our error bars if they show to be under-estimated. 
Please see SectionU!2]for details. 



3.1. MCMC 

We employed the MCMC method using the implementation de- 
scribed in [Gillon etal] ( |20 1 0| |20 1 2[ ) . In short, the radial veloci- 
ties are modeled using a Keplerian together with the prescription 
of the Rossiter-McLaughlin effect ( Rossiter||1924 McLaughlin 
| 1924[ ) provided by Gimenez1( |2006 l. The photometric model for 
eclipses (transits and occultations) is that of Mandel & Agol 
(2002), used without limb-darkening for occultations. The jump 
parameters are transit depth dF, impact parameter b, transit du- 
ration d, time of mid-transit Tq, period P, occultation depths 
dFocc for each wavelength and K2 - K^/l - e^p'/^ (where K 
and e denote the radial velocity semi-amplitude and eccentric- 
ity, respectively). The jump parameters Ve cos wand -s/esino) 
(where co denotes the argument of periastron) are used for the 
determination of the eccentricity. Limb darkening is accounted 
for by using the combinations ci = 2ui+U2 and 02 - ui - 2u2 of 
the calculated limb-darkening coefficients of [Claret & Bloenien 
( |2011| l following [Holman et aL]p006| l. With the exception of the 
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Fig. 2: The occultation observed with HAWK-I at 1.19 (im on 20 January 2010. The Hghtcurves obtained from the two jitter positions 
are shown. Top panel: the raw Hghtcurves together with the occultation and photometric model (left) and divided by the photometric 
model (right). Bottom panel; both hghtcurves corrected and normalized together with the best fit (left) and the residuals (right). 



a width equal to the error quoted by Claret & Bloemen (201 1 1) 



limb darkening parameters (for which we use a normal prior with et al. 2010[ l was fitted with the minimal model. Although we 

~ ~' ' '^^^^ derive photometric error-bars including scintillation, readout, 

background and photon noise, they are often under-estimated. 



we assume uniform prior distributions. We follow the method 
described by Enoch et al.| p010) using the mean stellar density, 
temperature and metallicity for determining the stellar mass and 
radius. In our whole analysis, we always ran at least two MCMC 
chains and checked convergence with the Gelman & Rubin test 
( [Gelman & Ru bin 1992). All time stamps are converted to the 
TDB time standard, as described by |Eastman et al.| ( i2010| . 



3.2. Photometric Model and Error Adaptation 

As described in Section [T] ground-based Hghtcurves are often 
affected by red noise correlated with external parameters. In 
our MCMC analysis, we have the possibility to include time, 
FWHM, coordinate shifts, and background variations in our 
model. This is done by multiplying the transit model by a poly- 
nomial (up to 4th degree) with respect to any combination of 
these parameters. The coefficients of the polynomial are not in- 
cluded as jump parameters in the MCMC but are found by min- 
imization of the residuals at each step. In order to account for 
airmass and stellar variability effects, we assume a second or- 
der polynomial with respect to time as the minimal accurate 
model for ground-based photometry. We checked more com- 
plex models by running MCMC chains of 10^ points on each 
lightcurve including higher orders of time dependence and addi- 
tional terms in FWHM, pixel position and background. Always, 
a more complex model was favoured over a simple one only if 
the Bayes factor (Schwarz||1978) estimated from the Bayesian 
Information Criterion indicated significantly higher probability 
(i.e. Z?i 2 > 100). The best photometric model functions are 
listed in Table [T] and are used in all subsequent analyses. For 
the H-band data, we kept the coordinate dependence described 
in |Anderson et al.| ( |2010| ), and the archive FTS lightcurve ( |Hebb| 



and we adapt them by accounting for additional white and red 
noise. For the white noise, we derive a scaling factor j6it from the 
ratio of the mean photometric error and the standard derivation 
of the photometric residuals. For the red noise, we obtain a scal- 
ing factor Pf by comparing the standard deviation of the binned 
photometric residuals to the standard deviation of the complete 



dataset, as described in detail by |Winn et aL ( 2008 1 and GiUon 



et al. (2010 1. Finally we multiply both scaUng factors to obtain 



the correction factors CF - li„ x for the photometric errors. 
In the subsequent analysis, all photometric error bars are multi- 
plied by these factors Anologously, we computed values for the 
radial velocity jitter which were added quadratically to the radial 
velocity errors. 



3.3. Summary of Tested Models 

Having derived the above factors, we analyzed the entire dataset. 
We did so by running chains of 10^ points on all photometric 
and radial velocity data. Next to the global analysis, we also per- 
formed analyses of subsets of hghtcurves (described in detail in 



Sections 4.1 and 



4.3.3 1. Additionally, we searched for any colour 



dependence in the transit depths by allowing for depth offsets 
between the different filters (Section [4. 3. 2| i and also derived in- 
dividual mid-transit times in order to check for transit timing 
variations (Section [4.3.2| i. In order to verify the result of |Burton] 
et al. (20121, we performed a global analysis in which we in- 
cluded their value, DFocc.z' - 880 +190 ppm as a Gaussian prior 
(Section |4.3.3[ ). The results are described in detail in Section 
|4] while all newly obtained hghtcurves are shown in Figure [T] 
(Transits), Figure |2](NB1 190 occultation) and Figure [3] (z' -band 
occultations). 
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Fig. 3: All occultation lightcurves obtained in the z' band. The upper six lightcurves were obtained with EulerCam and are unbinned, 
while the lower four lightcurves were obtained with TRAPPIST and binned in two minute intervals. Left panel: the raw lightcurves, 
together with the occultation and photometric model. Right panel: the lightcurves together with the occultation model, divided by 
the photometric model. The lower panel shows all data, corrected for the photometric model and binned in four minute intervals. 
Data and model are shown on the left, while the residuals are shown on the right side. 
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Fig. 5: The increase in photometric precision by combining up to 
nine lightcurves from TRAPPIST. The RMS in bin sizes of (from 
top to bottom) two, four, six, eight and 10 minutes are shown. 
The grey dashed lines show the expected 1 / -^n^ decrease for 
each time bin. 
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Fig. 6: The simultaneous transit observation performed with 
TRAPPIST (red) and EulerCam (blue). For clarity, the data are 
binned into 5 minute bins, and the models of the individual tran- 
sits are shown as solid lines. The binned data from the combi- 
nation of all 1+z' lightcurves (same as in Figure[8]l are shown in 
black for comparison. 



4. Results 

4.1. The TRAPPIST Transit Sequence 

In order to investigate the benefits of combining several 
lightcurves, we divided our set of nine 1h-z' TRAPPIST tran- 
sit lightcurves into subsets containing all possible combinations 
of one to nine lightcurves and performed an MCMC analysis on 
each of them using the procedure described in Section [3] We 
can observe how the solutions converge as we use an increasing 
number of transits in Figure [4] The outlier located at high tran- 
sit depth stems from the transit observed on February 10 2011, 
which is showing very high red noise. It is clearly visible that 
combinations favouring a large transit depth also require a larger 
transit width. This might be related to the presence of spots on 
the limb of the star. While they push the solution to lower transit 
depths, also slightly alter the measured time of ingress or egress. 
Figure|4]also shows the global solution and the solution obtained 
from the EulerCam transits alone. The results from TRAPPIST 
and EulerCam agree within their error bars, yet the EulerCam 
data find a slightly larger transit depth, width and impact param- 
eter. 

Next to investigating the parameters obtained from combi- 
nations of transits, one can also check the photometric improve- 
ment reached by combining an increasing number of lightcurves. 
In order to measure this effect, we use each set of combinations 
of n = 1 to n = 9 lightcurves, fold them on the best-fit period, 
and bin the data. Points before -0.05 and after -1-0.05 days from 
transit center are discarded in order to avoid phases that are not 
covered by all lightcurves. In Figure |5] we show the photometric 
RMS against the number of combined lightcurves and compare 
the decrease in RMS to the "best case" (decrease with 4=). The 

increase in precision is near that value, particularly if small time 
bins are used. By combining all nine TRAPPIST lightcurves, we 
obtain an RMS of 321 ppm for a moderately sized time bin of 5 
minutes. 



4.2. One Simultaneous Observation 

The transit of 23 January 201 1 was observed with EulerCam and 
TRAPPIST simultaneously, using an r'-Gunn filter on EulerCam 
and an I+z' filter on TRAPPIST. Figure [6] depicts the two 
lightcurves superimposed. It is obvious that the TRAPPIST data 
are showing an anomalously small transit depth and a short 
term brightening during the second half of the transit. With one 
lightcurve only, one might conclude that the planet crossed a 
star spot during transit. Since the EulerCam light curve was ob- 
served at shorter wavelengths and given the cooler temperature 
of the spot, we would expect the effect to be more pronounced 
here. However in the r'-band the feature is absent, excluding the 
spot-hypothesis. Removing the points during the second half of 
the transit in the TRAPPIST lightcurve, we obtain a value within 
1 cr of the EulerCam result. 

4.3. Global Analysis 

We performed a global MCMC analysis using all available radial 
velocity and photometric data including both, transits and occul- 
tations (described in detail in Section [3]). The results are shown 
in Table |2] 

4.3.1. Eccentricity 

[Anderson et al.|p010| l presented for WASP-19b a tentative non- 
zero eccentricity of 0.016;^q qq^ from the analysis of their H-band 
occultation. When including their data in our analysis, we de- 
rived a lower value for the eccentricity with a similar signif- 
icance 0.0077^gyy32- Rsmovlng the H-band data, we obtained 
0.006 l^gQQ43, clearly non-significant. We therefore do not see 
any evidence for a non-zero eccentricity of the WASP- 19 sys- 
tem. 

4.3.2. Transit Deptli and Timing Variations 

One of the checks we performed on the data was letting the 
transit depth vary for different filters, in order to search for 
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Fig. 4: The solutions found from fits to all possible subsets of TRAPPIST transits are shown. The single transits are shown in grey, 
combinations of three, five and seven transits are shown in orange, blue and green, respectively. The solution to all TRAPPIST and 
EulerCam lightcurves are shown in black and yellow, respectively. The global solution is shown in red. 



any large wavelength dependencies in the star/planet radii ratio 
which can be used to constrain models of the atmospheric trans- 
mission. The phase-folded and binned lightcurves for each filter 
are shown in Figure |8] and the respective radii ratios are shown 
in Figure|9]and listed in Table|3] The r', I-Hz' and z' band values 
match very well, only the IC value is slightly lower, 2.6 cr below 
the global solution. 

We also searched for any deviation from a linear ephemeris 
by performing a global analysis while fixing the ephemeris to 
the ephemeris derived from the global analysis but letting the in- 
dividual mid-transit times vary. The results are shown in Figure 
|7] and listed in Table |4] While there is some (expected) scatter 
around the linear ephemeris, none of the deviations exceed 2.7 cr 
and thus we find no evidence for TTVs in the WASP- 19 system. 



4.3.3. z'- band and 1.19 [xm Occultations 

We measure a z'- band occultation depth of 352 ±116 ppm from 
the combined analysis of the ten z'- band occultation lightcurves 
in our dataset. For individual lightcurves the occultation is well 
buried in the noise. In order to verify that our non-zero occul- 
tation depth is not caused by a systematic effect present in a 
single or a small number of lightcurves, we proceeded in a sim- 
ilar way as with the set of TRAPPIST transits (Section[4TT]i. We 
created all possible subsets containing at least five occultation 
lightcurves and analyzed them while fixing all parameters except 
the occultation depth to the values derived above. Histograms 
of the derived occultation depths are shown in Figure 10 The 



results obtained from fits of fewer light curves are consistent 
with the presented value although they have lower significance. 
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WASP- 19 



Jump parameters 


AF = (R„IR,f 


0.02018 ±0.00021 


b' = a* cos(i,,) [/?,] 


0.649 ±0.012 


Tu [d] 


0.06586t«™3 


r[0] - 2450000 [HJD] 


6029.59204 ±0.00013 


P[A] 


0.7888390 ± 2 x 10"' 


K2 [ms-'rf'/3] 


238.1 ±2.7 


^Focc-j [PPm] 


352 ± 116 


AF„„jvBii90 [ppm] 


1711!?« 
3216!ti 


AFo„// [ppm] 


cos (jj 


0.053 ± 0.020 


^[e sin 0) 


054+" °" 


Vv. sin/, cos p 


1 .85^A In 

-0. 1 y 


Vv, sin/, sin/? 


-0.27 ± 0.23 


Cl,r' 


1.123 ±0.040 


C2,r' 


-0.052 ± 0.037 


Cl.IC 


0.903 ± 0.034 


C2,IC 


-0.173 ±0.025 


C1.I+Z' 


0.840 ± 0.050 


C2,I+z' 


-0.251 ±0.053 


Cl.z' 


0.831 ±0.029 


C2.Z' 


-0.218 ±0.021 


Deduced parameters 


/f [ms-'] 


257.7 ± 2.9 




1.376 ±0.046 


Mp [MA 


1.165 ±0.068 


Pp [pA 


Q 447+0.027 

"•^^'-0.025 


e 


0077+° °'"'* 

-0.0032 


o) [deg] 




a [AU] 


0.01653 ±0.00046 


a/R. 


1.217 ±0.024 


ip [deg] 


79.54 ±0.33 




0.645 ± 0.012 




0.652 ±0.015 


r„„ - 2450000 [HJD] 


6030.77766 ± 0.00088 


yS [deg] 


-8 4+' ° 


r«, [K]" 


2058 ± 40 


logSp [cgs] 


3.184 ±0.015 


M. [Mq] 


968+° °'*'* 


^. IRq] 


0.994 ±0.031 


p. [p©] 


983+° °^'' 


"l.r' 


0.439 ± 0.020 


M2,r' 


0.246 ±0.015 


"l,I+z 


0.286 ± 0.026 


1'2,1+z 


0.268 ±0.019 


f l.IC 


0.326 ±0.016 


W2,IC 


0.2497 ± 0.0094 


"l,z' 


0.289 ±0.014 


W2.Z' 


0.2536 ± 0.0077 



Table 2: The median values together with the l-cr errors of the 
marginalized posterior PDF obtained from the MCMC analysis 
of all new data plus the radial velocities and lightcurves pub- 
lished b>J^bFeta^ and |Hellier 
|eFaL] ( |2011| i. 

Assuming an Albedo of A=0 and full redistribution from the 
planets day to night side, F=l. 



Filter 


Wavelength (width) [nm] 


Rp/R. 


r'-Gunn 


664 (99) 


0.1437 ±0.0013 


IC 


760 (91) 


0.1336 ±0.0026 


I+z' 


838 (190) 


0.14171 ±0.00094 


z'-Gunn 


912.4 (68) 


0.1428 ±0.0014 



Table 3: The planet/star radius ratios obtained from a global 
IVICMC analysis while allowing for filter-dependent transit 
depths. For the I+z' and z' filters, the widths given here have 
been determined from the combination of the filter transmission 
and detector response curves, for the r' and IC filters we give the 
equivalent widths. 



Epoch 


Mid-transit time [HJDtdb - 2450000] 


0-C [min] 


-1537 


4817.14633 ±0.00021 


-0.24 ±0.30 


-876 


5338.56927 ± 0.00023 


0.27 ±0.33 


-621 


5539.72327 ± 0.00030 


0.36 ± 0.44 


-583 


5569.69826 ± 0.00036 


-0.92 ± 0.53 


-564° 


5584.68693 ± 0.00024 


0.12 ±0.34 


-564* 


5584.68684 ±0.00019 


-0.00 ± 0.28 


-541 


5602.83138 ±0.00046 


1.78 ±0.67 


-536 


5606.77464 ± 0.00022 


0.44 ± 0.32 


-535 


5607.56241 ± 0.00033 


-1.10 ±0.47 


-516 


5622.55057 ± 0.00026 


-0.79 ±0.38 


-503 


5632.80612 ±0.00025 


0.14 ±0.36 


-474 


5655.68222 ± 0.00045 


-0.19 ±0.66 


-455 


5670.66976 ± 0.00064 


-0.78 ± 0.93 





6029.59250 ± 0.00035 


0.67 ±0.51 


43 


6063.51174 ±0.00030 


-0.54 ± 0.44 



Table 4: The mid-transit times and their deviation from a linear 
ephemeris. The values were obtained from the combined analy- 
sis of all transits while the ephemeris was fixed to the one quoted 
in Table |2] 

TRAPPIST Ughtcurve 

EuIerCam lightcurve 




-1500 



1000 -500 
Epoch 



Fig. 7: The O-C deviations of the individual transits from the 
ephemeris given in Table |2| are shown. The filled circle repre- 
sents the FTS transit of Hebb et al. (2010), the squares represent 
data obtained with EulerCam and the diamonds represent data 
from TRAPPIST 



Recently 'Burton et al.'('2012^ presented an occultation depth for 
WASP-19b of 880 ± 190 ppm based on one Ughtcurve obtained 
with ULTRACAM mounted at the NTT telescope at ESO La 
Silla observatory. We tried to reproduce this value by using it as 
a Gaussian prior (with a width equal the error bar on the mea- 
surement) in our MCMC. Even with this prior, the resulting oc- 
cultation depth is 466 + 97 ppm. We are thus not confirming the 



measurement of |Burton et al.| ( |2012| l but conclude that the occul- 
tation of WASP-I9b in z' band is significantly smaller 

In our global analysis we include the occultation observed 
with HAWK-I using the narrow-band NB1190 filter. From 
our data we measure the occultation depth at 1.19 |„im to be 
1711+^«Jppm. 
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Fig. 8: The combined lightcurves in each of the observed bands 
binned in five minute intervals. The lightcurves are a combina- 
tion of (from top to bottom) three, one, nine and two observa- 
tions. Using the displayed five minute intervals, the RMS of the 
residuals in the interval [-0.05,0.05] days are (from top to bot- 
tom) 432, 1021, 321 and 487 ppm. 
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Fig. 9: The planet/star radius ratios obtained from an MCMC 
analysis of all data while allowing for filter-dependent transit 
depths. The error-bars in wavelength are spanning the filters ef- 
fective width. Note that the deviation IC lightcurve is based on 
only a single observation. 



5. Discussion 

From the homogeneous set of TRAPPIST l+z' lightcurves, we 
can evaluate the photometric improvement obtained from com- 
bining lightcurves. As presented in Figure |5] we are not far from 
the ideal case of only white noise. Note that, even if combining 
as many as nine lightcurves, we are still gaining in photometric 
precision by adding additional transits. 

Following our observation strategy of deriving the most ac- 
curate measurement of the overall transit shape and thus the 
planetary parameters, we can reduce correlated noise, which 
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Fig. 10: A histogram of the results obtained from modeling com- 
binations of five (green), seven (blue) and nine (red) occulta- 
tions. 
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Fig. 1 1 : The model spectra of the dayside atmosphere of WASP- 
19b computed by 'Madhusudhan (2012 1 compared to observa- 
tions. The 



z'-band and 1.19 [xm observations presented in this 
work are shown in blue, while t he data of|Anderson et al.| ( |2010| 
|20n] ), [Burton et al.| ^20T2\ and'Gibs on et al.| ( |2010| l are shown 
in black. The two model atmospheres shown have been com- 
puted for the carbon-dominated (C/O = 1.1, red) and oxygen- 
dominated (C/0 - 0.4, green) case and are reproduced here with 
kind permission of N. Madhusudhan. 



can drastically affect single lightcurves. This is most evident 
if the same transit is observed simultaneously using different 
instruments. In the case of a combination of nine TRAPPIST 
lightcurves, we have measured the planet/star radius ratio with 
a precision of 0.7% in l+z' -band, while three observations 
with EulerCam in r' -band and a combination of two EulerCam 
lightcurves and a FTS yield slightly larger errors, giving a preci- 
sion of 0.9 and 1.0 %. While these values agree well within their 
error bars, a single lightcurve obtained in IC-band gives a 2.6 cr 
lower value. This lightcurve is showing a small flux increase dur- 
ing transit (possibly the signature of a star spot), and we suggest 
this possible radius variation to be verified with additional data. 
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Discarding this point, we see a flat optical transmission spectrum 
of WASP- 19b. Overall, the values derived from our analysis are 
in good agreement with the values which have been previously 
derived. 

With two new occupation measurements, we can proceed to 
constrain the chemical composition and structure of the plane- 
tary atmosphere. Madhusud han| ( [20T2 ) show models of WASP- 
1 9b for a carbon-dominated C/O- 1.1 and an oxygen-dominated 
C/O - 0.4 atmosphere. For oxygen-rich models, a strong absorp- 
tion around 0.9 ^m is expected from the TiO and VO leading 
to a smaller z'- band occupation depth. Following the occupa- 
tion measurement by Burton et al. ( 2012[ l,|Madhusudhan (2012') 
tentatively classified WASP- 19b as having a carbon-dominated 
atmosphere. Comparing our values to these models (see Figure 
11 1, we find our measurement of the z'- band matches extremely 
well the oxygen-rich model, showing a higher absorption, in- 
dicative of higher abundances in TiO and VO. These elements 
require a higher concentration of oxygen in the planetary atmo- 
sphere in order to contribute measurably to the planetary spec- 
trum. The 1.19 ^im value matches equally well the oxygen- and 
carbon-rich models. Thus, we suggest WASP-19b is a highly ir- 
radiated oxygen-dominated planet, fitting the 02 or upper 01 
Class defined by Madhusudhan (j2012j|. This should be con- 



firmed by a joint analysis of the previously published data and 
the measurements added in this work as the models of plane- 
tary emission are not unique. From the Spitzer data on WASP- 
19b ( [Anderson et al.||201 l| l, we know that WASP-19b does not 
show a temperature inversion. At first glance, this might seem 
unexpected for an oxygen-dominated atmosphere, as it is pre- 
cisely the molecules of which we measure higher concentrations 
that are causing temperature inversions. Still, as diff'erent wave- 
lengths probe different depths in the planetary atmosphere, the 
z'- band observations probe a deeper atmospheric layer, below 
the expected temperature inversion. As the star is known to be 
active, TiO and VO in the atmosphere might be destroyed by 
the elevated stellar UV flux, thus inhibiting a temperature inver- 
sion in the upper atmosphere, while lower in the atmosphere the 
intact TiO and VO are causing the measured absorption in the 
z' band. 



6. Conclusion 

We have carried out an in-depth observing campaign on WASP- 
19 collecting a total of 14 transit and 10 occupation lightcurves 
with EulerCam and TRAPPIST as well as one 1.19 ^m 
lightcurve with HAWK-1. From the large homogeneous set of 
nine TRAPPIST lightcurves, we demonstrate how both the at- 
tainable photometric precision and the accuracy of the derived 
parameters can be greatly improved by combining an increasing 
number of lightcurves. 

We have detected the z'-band occupation of WASP- 19b 
using Im class telescopes by the combined analysis of our 
lightcurves. We measure it at 352 + 116 ppm, more than a fac- 
tor of two smaller than previously published. From our HAWK-I 
data we obtain an occupation depth of 17 1 1^^730 ppm at 1.19 (,im. 
These results shed new light on the chemical composition of the 
planetary atmosphere, indicating a C/O ratio of C/O < 1, i.e. an 
oxygen-dominated atmosphere. 
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